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An overview is presented of a recently proposed "radically conservative" solu- 
Jl^ , tion to the problem of dark energy in cosmology. The proposal yields a model 

' universe which appears to be quantitatively viable, in terms of its fit to su- 

^ ' pernovae luminosity distances, the angular scale of the sound horizon in the 

cosmic microwave background (CMB) anisotropy spectrum, and the baryon 
acoustic oscillation scale. It may simultaneously resolve key anomalies relating 
to primordial lithium abundances, CMB ellipticity, the expansion age of the 
^ ' universe and the Hubble bubble feature. The model uses only general relativ- 

ity, and matter obeying the strong energy condition, but revisits operational 
00 ■ issues in interpreting average measurements in our presently inhomogeneous 

^\ • universe, from first principles. The present overview examines both the founda- 

ff^ ' tional issues concerning the definition of gravitational energy in a dynamically 

expanding space, the quantitative predictions of the new model and its best- 
fit cosmological parameters, and the prospects for an era of new observational 
tests in cosmology. 
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1. Introduction 

b . 

5^ , Dark energy is widely described as the biggest problem in cosmology today; 

one which may demand new physics and a theoretical paradigm shift. ^ In 
this paper I suggest that the solution to the problem of dark energy is inti- 
mately related to the correct understanding of observational anomalies - in 
particular, the observed abundance and emptiness of voids, which has also 
elicited separate calls for a paradigm shift. ^ I propose that the paradigm 
shift that is required to understand both these issues actually entails no 
"new" physics, but a revisiting of fundamental issues relating to the sub- 
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2 David L Wiltshire 

tlcty of the definition of gravitational energy within general relativity, and 
its relation to the careful modelling of the distribution of matter that we 
actually observe. 

The punchline is that cosmic acceleration can be understood as an ap- 
parent effect, and dark energy as a misidentification of those aspects of 
cosmological gravitational energy which by virtue of the strong equivalence 
principle cannot be localized,^ namely gradients in the quasilocal gravita- 
tional energy associated with spatial curvature gradients, and the kinetic 
energy of expansion, between bound systems and the volume-average posi- 
tion in freely expanding space. With this interpretation, a two-scale model 
can be constructed,^ and a simple exact solution'' of the relevant equations 
of cosmic evolution, the Buchert equations,^ can be found. This solution 
provides a fit to type la supernovae (Snela) luminosity distances which is 
statistically indistinguishable from that of the standard Lambda Cold Dark 
Matter (ACDM) model, while simultaneously satisfying key independent 
cosmological tests and offering the potential to resolve significant observa- 
tional anomalies.^ 

Over the past decade we have come to think of "dark energy" as a 
homogeneous isotropic fluid-like quantity, with a local pressure related to 
its energy density by P = wp, where w < —1/3, so that the strong energy 
condition is violated. In fact, observations appear to be most consistent 
with a pure vacuum energy, w = —I, and determination of the value of the 
parameter, w{z), as a function of redshift, z, is seen as the goal of "precision 
cosmology" . What is proposed here, if correct, will turn this situation on 
its head. We can look forward to an era of precision cosmology, but one in 
which the focus will be on the complex hierarchical structure of the universe, 
rather than any one simple fluid equation of state. Nonetheless, while the 
proposed solution is intimately connected to the growth of inhomogeneities, 
and their backreaction on the geometry of the universe,^ at its heart it 
addresses the question of the normalization of gravitational energy relative 
to observers within the observed structure. Thus I claim that the solution to 
the central foundational question does concern energy, and since "nothing" 
is "dark" the terminology "dark energy" is actually quite apt for the new 
proposal, if the community will allow the liberty of a change to the assumed 
definition of those words. 

In this paper I will give an overview of the proposed solution to the 
problem of dark energy, the extent of its present quantitative successes, 
and more importantly the directions for future work. I will present fewer 
technical details than may be found in papers already published^'^'^ or in 



February 2, 2008 15:23 



WSPC - Proceedings Trim Size: 9in x 6in dark 



Dark energy without dark energy 3 

preparation.*^''' I aim to give a general overview to researchers in astro- 
physics, particle physics or general relativity, who have no specific prior 
experience with the averaging of inhomogeneous cosmological models. 

2. The universe we observe 

Our most widely tested "concordance model" of the universe relies on the 
assumption of an isotropic homogeneous geometry at all epochs of cosmic 
evolution. By the evidence of the cosmic microwave background (CMB) 
radiation, the universe was very smooth at the time of last scattering, and 
these assumptions were completely well justified then. The departure from 
homogeneity was of order 6p/p^\0^'' in photons and the baryons that 
couple to them, and perhaps of order 5p/p~10~^ in non-baryonic dark 
matter, which gives rise to the potential wells responsible for the dominant 
Sachs- Wolfe effect. By the Copernican principle, the assumption of global 
isotropy and homogeneity is completely justified at the epoch of last scatter- 
ing, and it is safe to assume that the evolution of the universe was therefore 
extremely closely modelled by the Friedmann-Lemaitre- Robertson- Walker 
(FLRW) solutions at that epoch. Furthermore, if we consider the spectrum 
of CMB anisotropies^'^'^^ then overall it appears that globally the universe 
is very close to spatially flat. Its initial evolution at the time of last scat- 
tering would therefore have been very close to that of an Einstein-de Sitter 
model at that epoch, even in the case of the ACDM paradigm since "dark 
energy" only becomes appreciable at late epochs. 

At the present epoch, however, the distribution of matter is far from 
homogeneous on scales less than 150-300 Mpc. The actual universe has a 
sponge-like structure, dominated by huge voids. These voids are surrounded 
by bubble walls, and threaded by filaments, within which clusters of galaxies 
are located. Recent surveys suggest^^ that some 40-50% of the present 
volume of the universe is in voids of a characteristic scale 30/i~^ Mpc, where 
h is the dimensionless Hubble parameter, ffg = 100/ikm sec~^Mpc~^. If 
smaller minivoids^'^ and larger supervoids^^ are included, then our observed 
universe is presently void-dominated by volume. 

Quite apart from the fact that this observed structure appears emptier 
than the vistas that Newtonian iV-body simulations typically produce, the 
mere fact that the universe is presently inhomogeneous means that the 
assumptions implicit in the FLRW approximation can no longer be justified 
at the present epoch in the almost exact sense that they were justified at 
the epoch of last scattering. Homogeneity only applies at the present day 
in an average sense. The manner in which we take the average, and the 
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Fig. 1. Local voids and bubbles from the 6df survey. Courtesy of A. Fairall. 

operational issues associated with this are not trivial, since the problem of 
fitting the local geometry of bound systems to the full dynamical geometry 
of the evolving universe^^'^^ is a complicated one. 

Since a broadly isotropic Hubble flow is observed, an almost FLRW ge- 
ometry would appear to be a good approximation at some level of averaging, 
if our position is a typical one. I will demonstrate, however, that the as- 
sumption that the observed universe evolves exactly as a smooth featureless 
dust fluid means that we factor out central physically important questions 
which need to be understood to correctly relate our own observations to 
the average geometry. 



2.1. The Sandage de Vaucouleurs paradox 

There is a central foundational paradox concerning the expansion of the uni- 
verse, which others have called the "Hubble-de Vaucouleurs paradox" 
or the "Hubble-Sandage paradox" but which I will call the Sandage-de 
Vaucouleurs paradox since it was originally raised by Sandage and collabo- 
rators^" in objection to de Vaucouleurs' hierarchical cosmology.^^ 

The problem is that in the conventional way of thinking about cosmo- 
logical averages, below the scale of homogeneity we should expect large 
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statistical scatter in the peculiar velocities of galaxies. In fact, if we were 
to average on scales of order 20Mpc, which is about 10% of the scale of 
homogeneity, then the statistical scatter should be so large that no linear 
Hubble law can be extracted. Yet 20Mpc is the scale on which Hubble orig- 
inally extracted his famous linear law. By conventional understanding the 
statistical quietness of the local Hubble flow does not make sense. 

One might attempt to explain the Sandage-de Vaucouleurs paradox as 
a consequence of dark energy, since it is well-known that in any FLRW 
model which expands forever - with or without dark energy - peculiar ve- 
locities eventually decay. However, if one tries to numerically model 12.5 
Gyr of motion of very local galaxies, then it turns out that the initial con- 
ditions of each galaxy appear to have the most bearing on the problem.^'' 
In particular, one can fit a realistic motion with a cosmological constant, 
A, or alternatively in an open universe without A. Similar conclusions are 
also reached in constrained dark matter simulations of the Local Group. 
Furthermore, the ACDM parameters required for the velocity dispersion 
predicted by structure formation to match the observed velocity dispersion, 
do not coincide with the concordance parameters. Evidently a cosmolog- 
ical constant alone cannot explain the quietness and linearity of the local 
Hubble flow. 

A related, but not exactly equivalent, issue is that using a conventional 
kinematic approach, the peculiar velocities of local galaxies appear to be 
at least a few factors too small to have arisen from a smooth distribution 
of matter at the time of recombination.^^ 

3. Averaging, backreaction and dark energy: the debate 

The present distribution of matter is clearly very complex, and since we 
cannot solve the Einstein equations for this distribution of matter analyt- 
ically, there is an important question as to how we operationally match 
the average geometry of this distribution to the simple FLRW models that 

wc know how to solve. Given that the nonlinear growth of structure ap- 
pears to be roughly correlated to the epoch when cosmic acceleration is 
inferred to begin, a number of cosmologists have questioned whether the 
introduction of a smooth dark energy is a mistaken interpretation of the ob- 
servations.^ Attention has focused on the possibility that effects attributed 
to cosmic acceleration may actually be due to the backreaction from the 
growth of inhomogeneities in determining the geometry of the observed 
universe, without exotic dark energy. Different interpretations of a host of 
complex technical issues have led to a vigorous debate. 
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There are two large streams of research, which I wih not fociis on, for 
physical reasons. The first is backreaction in perturbation theory about 
an FLRW background, which became the focus of debate in 2005 following 
public attention generated by papers of Kolb and collaborators. ^''■^^ Within 
perturbation theory one may demonstrate that there is a potentially sig- 
nificant effect from backreaction.^'' However, this argument cannot be con- 
clusive - if the second-order terms in the perturbation expansion exhibit 
an effect which might be interpreted as cosmic acceleration, such an effect 
may go away when the third-order terms are considered, and so on. Pertur- 
bation theory is very relevant near the epoch of last scattering, when the 
assumption of homogeneity was extremely good; but by the present epoch 
the nonlinear structures are so numerous and complex that we are beyond 
the regime of its applicability. Thus perturbation theory cannot provide a 
complete solution, and will not be discussed here. 

The second stream of research I will not consider arc; those that involve 
exact inhomogeneous solutions of the Einstein equations, at the expense of 
introducing matter distributions which are unlikely. The spherically sym- 
metric Lemaitre-Tolman-Bondi (LTB) solutions are perhaps the most well- 
studied class of such models. While they may be very realistic descriptions 
of single voids, to apply them to the universe as a whole violates the Coper- 
nican principle, which I shall retain. One may obtain LTB solutions within 
which one can obtain reasonable fits to supernovae luminosity distances. 
However, in my view, given their high degree of symmetry, these are at 
best toy models, which one cannot hope to reproduce in structure forma- 
tion scenarios based on our knowledge of the power spectrum of density 
perturbations at the time of last scattering. 

To confront the actual inhomogeneous universe, which has no particu- 
lar spatial symmetries below the scale of homogeneity, we must deal with 
schemes that average the full non-linear Einstein equations. There are many 
schemes for constructing averages, including those of Zalaletdinov^* and 
Buchert.^ There are various grounds for favouring the approach of Zalalet- 
dinov,^^ which is fully covariant and averages all of the Einstein equations. 
However, Zalaletdinov's scheme is a general one, and for the cosmologi- 
cal problem additional assumptions are required. In Buchert's scheme one 
just average scalar quantities, and an additional integrability condition is 
required for the equations to close. However, with reasonable cosmologi- 
cal assumptions, the correlation tensor in Zalaletdinov's scheme takes the 
form of a spatial curvature,'^''' and Buchert's scheme can be realized as a 
consistent limit. Furthermore Buchert's scheme yields equations which 
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arc close to the Fricdmann equations. Given that the Friedmann equations 
have worked so well to date, Buchert's average would appear to give a nat- 
ural framework within which corrections to the Friedmann evolution can be 
quantitatively examined for the universe we observe. I shall adopt Buchert's 
scheme, with caveats to be discussed. 

Buchert's scheme strictly deals with irrotational dust cosmologies, char- 
acterized by an energy density, p{t, x), expansion, x), and shear, a{t, x), 
on a compact domain, I?, of a suitably defined spatial hypcrsurfacc of con- 
stant average time, and spatial 3-metric, ^Qij- Angle brackets are taken 
to denote the spatial volume averages, e.g., for the scalar curvature 

{T^)=(^j d3a;Vd^7^(^,x)^ /V{t) 

with V{t) = d^x-\/ det ^g. The important lesson of Buchert averaging is 
that time evolution and averaging to do not commute.^ Generally for any 
scalar 

^(*)-(^) = M-WW (1) 

The fact that the r.h.s. of (1) does not vanish, as is the case for the FLRW 
cosmologies, is a manifestation of backreaction. 

Applied to the equations of cosmic evolution one obtains the exact 

Buchert equations 

3^=87rG(p)-i(7e)-iQ, (2) 



a 



.a 



3- = -AnGip) + Q, (3) 



a 



St(p) +3-(p) =0, (4) 
where a{t) = [V{t)/V{t^f\ 

Q^l{{^^)-m-2{af (5) 
and following integrability condition follows from (2)-(5): 

dt{a^Q)+a''dt{a''{n))=Q. (6) 

It is observed that the kinematic backreaction term, Q, enters (3) with 

the same sign as a cosmological constant in the equivalent Raychaudhuri 
equation in the FLRW paradigm, but enters (2) with the opposite sign to 
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a cosmological constant in the Fricdmann equation. Given this fact, to- 
gether with the domain-dependence of the average, the question of whether 
Buchert's scheme can generate sufficient backreaction to give apparent ac- 
celeration for reaHstic initial conditions has been the subject of much de- 
bate. On reasonable grounds one might conclude"^^ that while the effects are 
real, they are too small in magnitude to give departures sufficiently large 
from the FLRW expectation to register as cosmic acceleration. 

It is at this point in the argument caution must be exercised. As every 
student of general relativity knows, one cannot simply write down a time 
parameter and assume that it is the parameter one measures on one's clock, 
without specifying how it is related to local invariants. We actually mea- 
sure luminosity distances^ and the deduction of cosmic acceleration involves 
two time derivatives. We therefore have to be extremely careful to opera- 
tionally specify how t is to be related to our own clocks. It must be observed 
that the scale factor a[t) is not related to an exact local metric, that has 
been substituted in Einstein's equations. Rather we must solve the Buchert 
equations, and determine how the best-fit almost-FLRW geometry that is 
obtained relates operationally to our own measurements. 

4. Finite infinity and gravitational energy 

In a completely arbitrary inhomogeneous universe the calibration of rods 
and clocks at one point relative to another can vary arbitrarily. However, our 
clocks and rods and those of stars in distant galaxies we observe, appear to 
a very good approximation to be determined by geodesies of ideal solutions 
- the Kerr and Schwarzschild geometries - in which space is asymptotically 
fiat, with an exact time symmetry at spatial infinity. This time symmetry, 
mathematically described by a timelike Killing vector, must be an idealiza- 
tion which breaks down at some level, since the universe is expanding. 

In his pioneering work on the fitting problem, EUis^^ suggested its so- 
lution should involve the notion of finite infinit-jf, ^^fi" , namely a timelike 
surface within which the dynamics of an isolated system such as the solar 
system can be treated without reference to the rest of the universe. After 
all, the matter in our galaxy and other typical galaxies broke away from 
the Hubble flow to be come a bound system well over 10 billion years ago. 
With suflacient computing resources, we can integrate the motions of stars 



^As finite infinity adds a new notion of infinity to the concepts timelike, spa- 
tial and null infinity, a new mathematical symbol, fi, is appropriate - in L^lijiX: 
\def\finfty{\mathop{\hbox{\it fi}}}. 
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within our galaxy for billions of years without considering the dynamics of 
the universe outside the galaxy. Thus within finite infinity spatial geometry 
might be considered to be effectively asymptotically flat, and governed by 
"almost" Killing vectors. 

The concept of finite infinity docs not appear to have been further devel- 
oped in the intervening two decades. However, given that the normalization 
of our clocks in the idealized Schwarzschild and Kerr geometries is related 
to the timelikc Killing vector at spatial infinity, finite infinity would appear 
to be the appropriate reference surface for the definition of gravitational 
energy, which in stationary spacetimes is tied to the asymptotically flat 
region.^^ In Newtonian terms, it is the scale at which we set the zero of the 
gravitational potential. 

The definition of gravitational energy in general is an extremely difficult 
problem, on account of the fact that space itself is dynamical and can carry 
energy and momentum. By the strong equivalence principle, since the laws 
of physics must coincide with those of special relativity at a point, it is 
only internal energy that can be localized in an energy-momentum tensor 
on the r.h.s. of the Einstein equations. Any uniquely relativistic aspects 
of gravitational energy associated with spatial curvature and geometrody- 
namics cannot be included in the energy momentum tensor, but are at best 
described by a quasilocal formulation.^^ 

Einstein himself worried about the problem of quasilocal gravitational 
energy, in terms of energy-momentum complexes, and many mathemati- 
cal relativists since Einstein have also considered the problem. It is quite 
possible that the general problem of a definition of quasilocal energy does 
not have a solution, since it depends on the split of space and time. Here 
I will not be interested in the general problem for an arbitrary solution 
of the Einstein equations, but in the specific problem for a universe which 
was effectively homogeneous and isotropic at last scattering, with nearly 
scale-invariant perturbations. 

The fact that quasilocal energy is not part and parcel of theoretical 
framework of the FLRW cosmology, is easily illustrated by the Newtonian 
perspective.^'* The l.h.s. of the standard Friedmann equation with A = 
can be regarded as the difference of a kinetic energy density per unit 
rest mass, E'kin = a^/(2a^) and a total energy density per unit rest mass 
Etot = —k/{2a?') of the opposite sign to the Gaussian curvature, k, while the 
energy-momentum tensor on the r.h.s. is the Newtonian potential energy 
per unit rest mass. The terms in the Einstein tensor represent forms of 
gravitational energy, but since they are identical for all observers in an 
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isotropic homogeneous geometry, one can always synchronize clocks and 
calibrate rods of ideal isotropic observers unambiguously. 

As soon as inhomogeneity enters the game, however, one will have gra- 
dients in both the kinetic energy of expansion, and in spatial curvature. 
Given the present value of the Hubble constant, it is likely that the spatial 
curvature gradient is the most significant effect. Quasilocal energy gradients 
between a bound system and finite infinity could conceivably be relevant for 
asymptotic galactic dynamics, or even the solution of the Pioneer anomaly. 
This speculation is left for future work. The proposal of Refs. [3,4] is con- 
cerned principally with dynamically evolving spatial curvature gradients. In 
this context, one should observe that the idea that negative spatial curva- 
ture is associated with positive gravitational energy, evident already in the 
Newtonian framework, remains true in the LTB models, where the Gaussian 
curvature is replaced by an inhomogeneous energy function E{r). 

Cosmological gravitational energy is largely uncharted territory in the 
more rigorous quasilocal framework^'^ within which the problem should ul- 
timately be framed. In any quasilocal formulation results depend crucially 
on the reference spacetime and surface of integration. Recently Chen, Liu 
and Nester'^"'' have obtained a result over which they expressed surprise, 
but which is consistent with the present proposal. They find that for an 
isotropic observer in synchronous gauge in a fc = — 1 Priedmann universe 
the quasilocal energy in their particular Hamiltonian formalism is nega- 
tive. A similar result is obtained using a difi'erent approach by Garecki.'^^ 
These results are expected in the current approach, since one is effectively 
subtracting a fiducial flat spacetime in each case, and the relative sign of 
energy depends on the observer. An isotropic fc = Friedmann observer 
has zero quasilocal energy in the approach of Chen, Liu and Nester; thus 
relative to the k = —1 geometry the k = Q geometry has negative quasilo- 
cal energy, but conversely relative to the A; = geometry the fc = — 1 has 
positive quasilocal energy. Our viewpoint here will be that the fiducial ref- 
erence point is the fc = geometry of the finite infinity region. This agrees 
with the Newtonian version of energy in the Friedmann equation, the LTB 
energy function, and with the idea that binding energy is negative. 



5. Finite infinity and the true critical density 

In the standard FLRW cosmology the critical density is proportional to the 
square of the Hubble parameter. However, in an inhomogeneous cosmology 
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this is no longer true, 

^ 8^ 

as the average Hubble parameter, ffav, is related to the spatial scale of 
averaging and the clock chosen. In the spirit of the spherical collapse model, 
we are familiar with the idea that a perturbation near the time of last 
scattering can evolve as an FLRW model of different density, until it turns 
around and collapses, by which point other regions closer to the true mean 
density are within the past light cone. 

In a universe which grew from a nearly scale-invariant spectrum of 
perturbations, there will always be spatial correlations of some density per- 
turbation which is relatively far from the mean, once one samples on scales 
close to the spatial scale of the particle horizon. This is a simple conse- 
quence of cosmic variance, and is illustrated by the following analogy. Take 
a bucket filled with grains of sand of identical diameter, whose density is 
Gaussian distributed about a mean. Provided the container is much larger 
than the individual grains then the mean density of the bucket can be ex- 
pected to be extremely close to the mean density of all the sand from the 
beach on which it was collected. If one repeats the exercise for stones of 
larger and larger diameter, the mean density of the bucket will on average 
differ more and more from the mean density of the beach as the diameter of 
the stones becomes comparable to the size of the bucket, by a a/TV statistic. 

Cosmic variance'' in the actual universe is complicated by three things. 
Firstly, since general relativity is causal, the bucket in the analogy is the 
past light cone which grows with time. Secondly, the perturbations are not 
uniform density but contain other perturbations embedded within them, 
like Russian dolls. Thirdly, below a scale of apparent homogeneity the initial 
perturbations can no longer be thought of as perturbations. They have 
undergone a non-linear transition to form the structure that we see: stars, 
globular clusters, and clusters of galaxies strung in filaments and bubbles 
around voids. 

Every inflationary cosmologist who has thought about the foundations 
of the subject realises that the present density of the observable universe 
need not be the same as the whole universe, observable and unobservable, 
and in general these should differ on account of cosmic variance in a scale- 



''I use the term cosmic variance in a general sense, rather than restricting it to one 
statistical aspect of the CMB temperature fluctuations which arise as an observable 
consequence of the variance in the density perturbation spectrum. 
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free density perturbation spectrum. Nonetheless, in the FLRW paradigm 
we unwittingly end up violating this. As long as we assume that the average 
density evolves by the Friedmann equation, we implicitly assume that the 
mean density in our past light cone today is the same as the whole ensemble 
of observable universes. As long as we evolve the average geometry purely 
by the Friedmann equation, with equality in (7), we will miss an important 
aspect of actual cosmic evolution. Furthermore, as long we study inhomo- 
geneity in the context of a Swiss cheese modcl^^ by cutting and pasting holes 
into a background which still evolves by the Friedmann equation alone, we 
will still fail to uncover essential features of cosmic evolution. 

There is perhaps no obvious way to define finite infinity in an arbitrary 
inhomogeneous background. To proceed I make the crucial observation that 
since our universe was effectively homogeneous and isotropic at last scat- 
tering, a notion of a universal critical density scale did exist then. It was 
the density required for gravity to overcome the initial uniform expansion 
velocity of the dust fluid. I will assume, as is consistent with primordial 
inflation, that the present horizon volume of the universe was very close to 
the critical density at last scattering, with scale-invariant perturbations. 



In view of the existence of backreaction, what is required is a notion of 

true critical density without assuming evolution via the Friedmann equa- 
tion. This should be a dynamically evolving demarcation scale between sys- 
tems which will become bound and those that are unbound, given available 
data within the past light cone at any epoch. As outlined in Ref. [3], and 
depicted in Fig. 2, finite infinity is defined in terms of a scale over which the 
average expansion, including virialized, collapsing and marginally expand- 
ing regions is zero, (t?)^ = 0, while being positive outside. Then p„ = 
Finite infinity is a non-static boundary analogous to the spheres cut out 
in the Swiss cheese model, but it involves average geometry rather than 




e>o 



Finite infinity <0>=O 



Fig. 2. Finite infinity, fi. 
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matching exact solutions, and no assumptions about homogeneity are made 
outside finite infinity. Finite infinity represents a physical scale expected to 
lie outside virialized galaxy clusters, but within the filamentary walls sur- 
rounding voids. Since it is a scale related to the true critical density, space 
at finite infinity boundaries can be described by the spatially flat metric 

dsj = -dr^ + aj (tw) [drjl + r?^ dfJ^] , (8) 

where df2^ = dO^ + sin^ 9 d(l)^ . While I use the term "bubble wall" it should 

be recognized that an ideal bubble wall would consist of finite infinity re- 
gions touching exactly at their boundaries. Due to the existence of minivoids 
one should take care in identifying the scale of local walls empirically. 



6. Average homogeneity in a lumpy universe 

An important issue that has to be faced is the split of space and time. Den- 
sity is not a covariant quantity, but will depend on the foliation of hyper- 
surfaces chosen to average on, and the specification of the observers within 
the hypersurfaces. The hypersurfaces chosen by Buchert are the standard 
ones based on surfaces of constant, t, where this is the affine parameter on 
geodesies of ideal comoving observers. This choice, which is analogous to 
the choice of synchronous gauge in FLRW models, is always possible for an 
energy-momentum tensor with dust in the absence of vorticity. 

While it is entirely reasonable that Buchert 's choice of gauge is consis- 
tent at the volume-average position, provided one averages over suitably 
large spatial volumes, there is a problem of physical interpretation since 
observers in bound systems are located in places where the physical as- 
sumptions in Buchert 's gauge choice no longer apply. As soon as regions 
start to collapse, geodesies will cross, vorticity comes into play and comov- 
ing coordinates cannot be chosen globally. Actual observers are in locations 
where the geometry is better modelled by the vacuum Schwarzschild and 
Kerr solutions, rather than by a rotationless dust fluid. 

Since a nearly isotropic Hubble flow is observed, it is clear that at some 
level a uniform Hubble flow should be obtained, despite the large scale 
inhomogeneity. I will make a choice of surfaces of homogeneity which im- 
plicitly solves the Sandage-de Vaucouleurs paradox by the assumption that 
on suitably small scales of averaging the bare Hubble flow is uniform below 
the scale of apparent homogeneity, 

' ^^"^'^-'-W^ =l{e)^ =--- = m, (9) 



£r{t) dt 3' '^1 3 
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t being the local proper time in each averaging region Vj, and ir = V^/^ = 

1/3 

[Jp d"^3;-\/^] a relevant proper distance. Here no Dj can be contained 
within a finite infinity region. The important point is that although it will 
appear to any one observer that voids expand faster than bubble walls, the 
clocks of observers within voids will also tick faster than those of observers 
at finite infinity on account of the relative gravitational energy gradient. 
Thus, provided the gravitational energy and spatial curvature gradients 
are correlated, apparent variance in the Hubble flow for any one observer is 
nonetheless perfectly consistent with uniformity of the bare, or "quasilocally 
measurable" Hubble flow. 

The uniform expansion gauge, which is equivalent to one of the stan- 
dard gauges in cosmological perturbation theory,^^ is given as an ansatz for 
defining the surfaces of homogeneity. However, it should be noted that this 
ansatz is less restrictive than the standard ansatz implicit in the FLRW 
models, where expansion is uniform, and all ideal isotropic observers also 
have synchronous clocks and measure the same spatial curvature on the rel- 
evant surfaces of homogeneity. Here uniformity in the bare Hubble expan- 
sion is retained, but spatial curvature and gravitational energy of isotropic 
observers vary in a correlated fashion. 

It should be emphasized that the Copernican principle is retained here. 
If we compare ideal isotropic observers, some in galaxies and some in voids, 
they will also measure an isotropic CMB. However, they will potentially 
measure different mean CMB temperatures, and different angular scales in 
the CMB anisotropics. These two observations respectively relate to relative 
gravitational energy and relative spatial curvature. While the clocks of all 
ideal observers are synchronized initially, when the FLRW approximation 
was a good one, ultimately they will diverge; by 38% on average at the 
present epoch. Since we exchange photons with other bound systems, which 
are within finite infinity regions, such a large clock-rate variance is not 
directly detectable in experiments conceived to date. 

The uniform bare Hubble constant gauge is more than simply an ansatz 
which implicitly solves the Sandage-de Vaucoulers paradox while dealing 
with the question of cosmological gravitational energy. As is argued in a 
separate paper,* it may be realized as a consequence of the strong equiva- 
lence principle applied to expanding space. Among the class of all possible 
motions of a timelike geodesic congruence there is a class of conformal mo- 
tions which do not isolate any direction in space as preferred, in contrast to 
individual boosts. The statement then is that for such motions we cannot 
distinguish the circumstance in which the particles are moving from a com- 



February 2, 2008 15:23 



WSPC - Proceedings Trim Size: 9in x 6in dark 



Dark energy without dark energy 15 



mon origin in a static Minkowski space, or alternatively are at rest in the 
underlying expanding universe. The two situations are equivalent; giving 
a Cosmological Principle of Equivalence.^ In terms of the conceptual basis 
of relativistic cosmology, this amounts to a refinement of the Principle of 
Inertia. Not only can we not distinguish which one of Galileo's ships'^, or of 
Einstein's trains is the one that is moving; for conformal motions we can 
also not distinguish whether it is the collection of ships or trains that are 
moving, or the 'sea' or 'railyard' expanding in between. 

7. The fractal bubble model 

In Ref. [3] I have written down a Buchert average of the Einstein equations 
based on two scales, finite infinity regions within which the geometry is 
given by (8), and voids within which the geometry is negatively curved 
with local scale factor, The geometry near the centres of voids is given 



The Buchert equations are not written in terms of the local geometry (8) 
or (10) at cither finite infinity or the void centres, but in terms of an inter- 
mediate volume-average location with a volume-average time parameter, 
t, and an average scale factor, a, defined by 



/vi and /wi = 1 — /vi being the respective initial void and wall volume 
fractions at last scattering. Initially /vi 4C 1, and /wi — 1, since the universe 
is homogeneous and isotropic at last-scattering, evolving like an Einstein 
de Sitter one. The Buchert average is constructed over the entire present 
particle horizon volume. 

Although the actual surfaces of homogeneity will not coincide with 
Buchcrt's ones on small scales, we can still make use of Buchcrt's scheme, 
provided that we take care in defining the relationship between the volume- 
average quantities of Buchert's scheme, and those in a finite infinity region. 
Ultimately, in order to deal with actual gradients in spatial curvature be- 
tween finite infinity and void centres, it may be necessary to use a scheme 
such as that of Zalaletdinov,^* or to look at the problem in terms of Ricci 



by 



d4 = -dT^ + Ov'(rv) [d7?2 + sinh2(ryv)dfi2] . 



(10) 




(11) 



'^Actually Galileo just talked about a single ship, in a context more like Einstein's closed 
elevator, but the principle is the same. 



February 2, 2008 15:23 



WSPC - Proceedings Trim Size: 9in x 6in dark 



16 David L Wiltshire 

flow.'''' This may be necessary, for example, in determining the most ap- 
propriate gauge for the background to structure formation simulations. The 
approach adopted in Refs. [3,4,6] is that the Buchert equations are solved in 
volume average time, but the uniform expansion ansatz is used to calibrate 
observable quantities relative to observers at finite infinity. This involves 
a dressing of cosmological parameters over and above the volume factors 
considered by Buchert and Carfora.^^'^'^ 

The independent Buchert equations, including the integrability condi- 
tion that ensures their closure, may be written 

^M + ^fc + ^e = i' (12) 
a-^dt (p.^R^a^^ + a-^dt (fl^^H^a^'^ = , (13) 

where the volume-average matter, curvature and kinematic backreaction 
parameters are respectively 

8ttGp -k f -"^/^f 1/3 _ / 2 
o _ z__LM0Z0 ■ o - ^v-'v' o /v (TA\ 

The average curvature is due to the voids only, which are assumed to have 
fcv < 0, an overdot denotes a derivative w.r.t. volume-average time, t, and 
H = a/a is the volume-average or bare Hubble parameter. 

Although the quantity, H, has no particular significance for a general 
Buchert average, for our particular assumptions it represents the underlying 
uniform quasilocally measured Hubble flow. The quantities Cl^^, fi^ and 
CIq are then interpreted as the bare cosmological parameters, relevant to 
a comoving isotropic observer at an average position in freely expanding 
space, which will be within a void, but not at its centre. These quantities 
are the closest analogues of the standard FLRW density parameters. 

It must be recalled that the scale factor a does not correspond to an 
exact geometry substituted into the Einstein equations, and then solved. 
Rather we integrate the Buchert equations and then reconstruct the average 
spherically symmetric geometry 

ds^ = -df + a^{t) df + Ai-n, t) dn^, (15) 

the area function A being defined by a horizon-volume average.^ The fact 
that (15) is spherically symmetric reflects the fact that the geometry is re- 
constructed by taking an average on our radial null geodesies. It is therefore 
not an LTB model, since exact spherical symmetry has not been assumed 
in the Einstein equations. 
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Since (15) does not correspond to the local geometry of a finite in- 
finity observer, if we assume tfiat our rods and clocks differ little in cal- 
ibration from those at finite infinity, we still have to be careful in re- 
lating our local geometry (8) to (15). We do this by conformal match- 
ing of the radial null geodesies of (8) and (15), which requires that 
d'yw = /wi^''''d?7/[7 (1 — /v)^^"^]. It then turns out that the finite infinity 
geometry may be rewritten 

dsj = -drl + 32 [df + rlifj, Tw) dn^] (16) 

where = 7(1 — fv)^^"^ fwi^'^^^Vwifj^T^v)- The geometry (16), which now 
again has a more general spherically symmetric form, is effectively the clos- 
est match to the FLRW geometry we usually attempt to fit to the whole 
universe with the assumption that spatial curvature everywhere matches 
our own, and clocks everywhere arc synchronized to our own wall time, t^. 

In place of the bare cosmological parameters (14), one can define con- 
ventional dressed parameters with respect to the geometry (16), relevant 
to "wall observers" such as ourselves. In particular, the dressed density 
parameter, Qj^, defined according to = 'y^^j^^ is effectively the con- 
ventional parameter, whose numerical value will be similar to that which 
we infer in FLRW models. Since (16) is not an FLRW metric it does not 
make particular sense to supplement by additional dressed parameters, 
as they will not sum to unity. 

What is more significant cosmologically is that the dressed geometry 
(16) will yield a dressed Hubble parameter, 

Ido 
dr™ 



^=-3^=7^-3*7 = 7^-7-'! 7- (17) 



where a = a/7, which differs from the bare Hubble parameter, H. Sim- 
ilarly, a wall observer will determine a dressed deceleration parameter 
q = —d/{H^a), which differs from the bare deceleration parameter q = 

-'ti/{ifa). 

The general solution of the two scale Buchert equations was found in 
Ref. [4]. The general solution is specified by four independent parameters, 
i/p, f^^g, /vo and Ei, which are respectively the bare Hubble constant, the 
bare matter density, the present epoch void volume fraction, and a small 
parameter ej <C 1 related to the initial void fraction. However, two of the 
four parameters are greatly restricted by taking priors^ at the surface of 
last scattering consistent with the CMB. Furthermore, there is a tracker 
solution to which all solutions with these priors approach to within 1% by 
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a redshift z ~ 37. This effectively reduces the number of free parameters to 
two; we take these to be and /vo- The tracker solution is given by 

3/vo5oi+(l-/vo)(2 + /vo) (18) 



2 + /vo 



f = ^-^""-^o^ (19) 

3/v0^o*+(l-/vo)(2 + /vo) ' 

in terms of the volume-average cosmic time, t. Since the lapse function and 
bare Hubble parameter are related to the void fraction by 7 = 1 + i/v = 
^Ht, the relation drw = 7^^d< may be readily integrated to obtain the wall 
time, Tw, relevant to observers in galaxies. 



2(l-Ao)(2 + /.o) / 9UHot \ 

^-'^ 21UHo ^"l,'+2(l-Ao)(2 + Ao)J- 

As is observed in Ref . [4] , expressions for many relevant observable quan- 
tities, including effective dressed luminosity and angular diameter distances, 
dj^ and d^, are readily obtained. It is particularly interesting to compare 
the bare volume-average tracker solution deceleration parameter. 



2 (1 - /v 



i2 



(21) 



" (2 + /v)2 ' 
with the corresponding dressed deceleration parameter 

„ _ - (1 - /v) (8/v^ + 39/v^ - 12/v - 8) 
{A + U + AU^) 

At early times, when /v <C 1, both deceleration parameters begin at the 
Einstein-de Sitter value, q q i. The bare deceleration parameter al- 
ways remains positive, meaning that a volume average observer in freely 
expanding space detects no cosmic acceleration, in accord with the intuition 
of the critics. However, the dressed deceleration parameter changes sign 
at epoch when ~ 0.5867, at a zero of the cubic in (22), so that a wall 
observer does detect apparent acceleration. 

The fractal bubble (FB) model provides a solution to the observational 
coincidence that the onset of cosmic acceleration appears to roughly coin- 
cide with the growth of large structures. Cosmic acceleration is an apparent 
effect that arises from gravitational energy gradients related to spatial cur- 
vature gradients. It begins when the void fraction reaches 59% at a redshift 
z ~ 0.9. This statement itself should be able to be tested in future. Tradi- 
tionally observational cosmologists define voids in terms of a density con- 
trast threshold, when quoting void volume fractions. For example, one 
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might argue as to whether one should take Sp/p < —0.9 or some other 
bound. In the present model, the void volume fraction is related technically 
to those regions within the past light cone which are not within finite in- 
finity. With some further refinement, it should be possible to relate this to 
an observational definition of void fractions. Then as surveys improve, the 
claims about the properties of /v(^) can be compared with observation. 

This illustrates the power of the FB model as compared to the spatially 
flat ACDM model. Both models depend on the Hubble parameter and one 
other free parameter. However, is not directly observable, whereas the 
void-volume fraction, /v, is empirically observable in principle. 

8. Observational tests 

In Ref. [6] we have examined the fit of the FB model to the RiessOT gold data 
set^^ of Snela. We find that with 182 data points and two degrees of freedom 
the best-fit = 162.7, i.e., a of approximately 0.9 per degree of freedom, 
which is a good fit. While a marginally lower = 158.7 may be obtained for 
the best-fit flat ACDM model (with = 62.6km scc-iMpc~\ Qj^^^ = 
0.34, = 0.66), the difference is not significant. In fact, a Bayesian 
model comparison of the FB model against a flat ACDM model with priors 
55 < iJp < 75km sec"'^ Mpc"'^, 0.01 < i^j(^Q < 0.5. gives a Baycs factor of 
\nB = 0.27 in favour of the FB model. Statistically, such a small margin 
is "not worth more than a bare mention" or "inconclusive";^^ that is to 
say, the two models are statistically indistinguishable. 

The residual difference A/i = /ipg — /^cmpty, m the standard dis- 
tance modulus, p = 51ogiQ(d^) +25, of the best-fit FB model from 
that of a coasting Milne universe of the same Hubble constant, ~ 
61.7 km sec~^ Mpc~^, is plotted in Fig. 3 and compared with binned data 
from the Riess07 gold data set. Apparent acceleration occurs for positive 
residuals in the range, z < 0.9. (Note, however, that the exact range of red- 
shifts corresponding to apparent acceleration also depends on the value of 
the Hubble constant of the Milne universe used to compute the residual.) 
The equivalent theoretical residual for the spatially flat ACDM model with 
the same values of and fi^Q is also shown. The difference in the gra- 
dients of the two curves should lead to observable differences if observing 
programs which are trying to measure the "equation of state, P = wp, 
w{z), of dark energy" can achieve sufficient accuracy. 

In Fig. 4 we address three major independent cosmological tests. In ad- 
dition to statistical confidence limits in the (^q, ^^mq) para-mctcr space we 
also consider the angular scale of the sound horizon in the CMB anisotropy 
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Fig. 3. The difference in the distance modulus, = 51ogj()((i^) + 25, with in units 
Mpc, of the FB model with //(j = 61.7km sec^-"- Mpc^"', ^j^jq = 0.326 from that of an 
empty coasting Milne universe, with the same value of //jj. The solid curve shows the FB 
model expectation, and the dashed curve the expectation for a spatially flat ACDM model 
with the same values of H^^ and f^^^y Whiskers indicate how the statistical uncertainties, 
shown as boxes, move when the background value of H^^ for the Milne universe which is 
subtracted is varied within the 2(t limits. For further details see Ref. [6]. 

spectrum, ^"^'^^ and the comoving scale of the baryon acoustic oscillation 
(BAO) seen in galaxy clustering*^ using the relevant dressed geometry (16). 

Ideally we should recompute the spectrum of Dopplcr peaks for the 
FB model. However, this is a mammoth task, as the standard numerical 
codes have been written solely for FLRW models, and every step has to be 
carefully reconsidered. For this reason we first ask whether parameters exist 
for which the effective angular diameter scale of the sound horizon matches 
the angular scale of the sound horizon, 5 = 0.01 rad, of the ACDM model, 
as determined by WMAP.^° Since there is no change to the physics of 
recombination, but just an overall change to the calibration of cosmological 
parameters, this is entirely reasonable. In Fig. 4 we plot parameter ranges 
which match the 5 = 0.01 rad sound horizon scale to within 2%, 4% and 
6%, using the calculation of the sound horizon given in Ref. [3], Sec. 7.2. 
The 2% contour would roughly correspond to the 2(t limit if the WMAP 
uncertainties for the ACDM model are maintained. As this can only be 
confirmed by detailed computation of the Doppler peaks, the additional 
levels have been chosen cautiously. 

In the case of the BAO scale, as we do not yet have the resources to 
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Fig. 4. Various cosmological tests in the parameter space of the dressed Hubble con- 
stant, Hg, and the dressed matter density parameter, ^j^^q- (a) Statistical confidence 
limits on Snela in the Riess07 gold data set;*^ (b) Parameters which match S = 0.01 
angular scale of sound horizon^^'^^ of ACDM model to within 2%, 4%, 6%; (c) Param- 
eters which match effective comoving 104^"^ Mpc scale of baryon acoustic oscillation** 
to within 2%, 4%, 6%; (d) Overlay of panels (a),(b),(c). 



analyse the galaxy chistering data directly, we also begin with a simple 
but effective check. Since the dressed geometry (16) does provide an ef- 
fective almost-FLRW metric adapted to our clocks and rods in spatially 
flat regions, the effective comoving scale in this geometry should match the 
corresponding observed BAO scale of 104ft,~'^Mpc. In Fig. 4 we therefore 
also plot parameter values which match this scale to within 2%, 4% or 6%. 

It should be noted that even if Snela are disregarded, the param- 
eters which fit the two independent tests relating to the sound hori- 
zon and the BAO scale agree with each other, to the accuracy shown, 
for values of the Hubble constant which include the value of Sandage 
et al.f^ but not for the values of greater than 70 km sec~^ Mpc~ 
which best-fit the WMAP data^^-" with the FLRW model. The value 

= 62.3 ±1.3 (stat) ± 5.0 (syst) km sec~^ Mpc~^ determined by the Hub- 
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ble Key Team of Sandage et a/.^^ has been controversial, given the 14% 
difference from values which best-fit the WMAP data with the ACDM 
model. ^"'^^ However, the WMAP analysis only constitutes a direct mea- 
surement of CMB temperature anisotropics; the determination of cosmo- 
logical parameters involves model assumptions. In the FB model, the model 
assumptions are different. In particular, on account of the differences in 
gravitational energy and local spatial curvature between observers in bound 
systems, and those at the volume average in freely expanding space, the cal- 
ibration of quantities related to the CMB has to be carefully reconsidered. 
Observers at the volume average detect a cooler mean CMB temperature^ 
Tq = 7~^Tq = 1.98 K than the Tq = 2.73 K we measure. Accounting for 
these differences leads to concordance for parameter values which include 
the Hubble constant of Sandage et al.^^ 

Table 1. Best-fit cosmological parameters derived from the independent parameters, 

dressed Hubble constant 
present void volume fraction 
mean lapse function 
bare density parameter 
conveiitioual dressed density parameter 
mass ratio of 

non-baryonic dark matter to baryonic matter 
bare Hubble constant 
effective dressed deceleration parameter 
age of universe measured in a galaxy 

Note: 1 (7 statistical uncertainties from Snela alone are shown. 



= 61.7Ii 1 km sec-i Mpc"! 

/vo = o.76iEj:J2 

;^ _ 1 001 +0.061 
''mo ~ "-^^^3-0.069 

o — n -j-j+o ii 

(pMO ^Bp)/^SO ~ '^•^^2.4 

= 48.2I2 ^ km sec-1 Mpc'^ 
= -0.0428l°;°J?5 
Two = 14.7t°i Gyr. 



Best-fit cosmological parameters are shown in Table 1. As statistical 
bounds are not available for the sound horizon and BAO tests, we show la 
uncertainties from Snela only. However, given that there are many inde- 
pendent estimates of the dressed matter density, ^^j^/o' ^^'^ expect that the 
uncertainties quoted can be significantly reduced if such priors are imposed. 

8.1. Resolving the lithium abundance anomaly 

The angular scale of the sound horizon and the BAO tests have been applied 
assuming a volume-average baryon-to-photon ratio in the range ry^^ = 4.6- 

5.6 X 10"^° adopted by Tytler et al^^ prior to the release of WMAPl. With 
this range it is possible to achieve concordance with lithium abundances. 
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while also better fitting helium abundances. This may resolve the primordial 
lithium abundance anomaly. ^^'^"^ 

With the 2003 WMAPl release/*^ the baryon-to-photon ratio was in- 
creased to the very upper range of values that had previously been consid- 
ered, due to its effect on the ratio of the heights of the first two Doppler 
peaks. This ratio of peak heights is sensitive to the mass ratio of baryons to 
non-baryonic dark matter - rather than directly to the baryon-to-photon 
ratio - as it depends physically on baryon drag in the primordial plasma. 
The fit to the Doppler peaks required more baryons than the range of Tytler 
et al^^ admitted, when calibrated with the FLRW model. In the FB cali- 
bration, on account of the difference between the bare and dressed density 
parameters, a bare value of fi^j, ~ 0.03 nonetheless corresponds to a con- 
ventional dressed value VL^^ ~ 0.08, and an overall mass ratio of baryonic 
matter to non-baryonic dark matter typically of about 1:3, which is larger 
than for ACDM. This would certainly indicate sufficient baryon drag to 
accommodate the ratio of the first two peak heights. 

8.2. Spatial curvature and the ellipticity anomaly 

Since the release of the Boomerang experiments^ in 2000, it has been gener- 
ally assumed that the angular position of the first Doppler peak is a measure 
of the spatial curvature of the universe. However, this inference assumes 
that the spatial curvature is the same everywhere, as is appropriate for the 
FLRW models. In the present model there are spatial curvature gradients, 
and we must revisit the calculation from first principles, as outlined in Ref. 
[3]. Insofar as the angular scale of the sound horizon is reproduced for the 
FB calibration, the overall angular scale cannot be regarded as a measure 
of average spatial curvature. 

If there is a sizable negative average spatial curvature at the present 
epoch, then there must be ways of detecting it other than via the mea- 
surement of the angular scale of the Doppler peaks. Such effects do indeed 
arise when one considers more subtle mcasiirements associated with the 
average geodesic deviation of null geodesies. Indeed one prediction is that 
there should be nontrivial ellipticity in the CMB anisotropics on account of 
greater geodesic mixing. This effect is in fact observed,''^ and is an impor- 
tant anomaly for the standard ACDM paradigm, overlooked by the majority 
of the community to date. 

A detailed quantification of the degree of ellipticity expected, and its 
comparison with the observed signal, will be an important test of the FB 
model in future. 
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8.3. The WMAP3 normalization 

It has been recently noted that the vahies of the normahzation of the 
primordial spectrum (7g~0.76 and matter content rij^^Q~0.24 implied by 
WMAP3 are barely compatible with the abundances of massive clusters 
determined from X-ray measurements.^"^ In fact, Snela also best fit the 
ACDM model with higher values of ^j^^ comparable to the the best fit 
value for the FB model, namely ^j^q ~ 0.33. In the FB model the best-fit 
void fraction, /vo ^ 0.76, appears to be the quantity that is mimicked by 
dark energy fraction, fi^^ , as far as the WMAP normalization to the ACDM 
model is concerned. The flat ACDM model constrains i^^g = 1 — f^^o' ^'^^ 
the FB model, the corresponding constraint, ^j^q ~ i(l — /vo)(2 + /vo)) 
gives quite different larger values, consistent with many other observational 
determinations of the conventional matter density parameter. 



8.4. The expansion age 

Structure formation scenarios in standard ACDM model have some diffi- 
culty in explaining the observed apparent very early formation of galaxies. 
Of course, structure formation contains many model dependent assump- 
tions, so direct measurements of the ages of metal poor stars in old globu- 
lar clusters are ultimately a better indicator. The observational bounds are 
not very tight, because of the large uncertainties involved in particular nu- 
clear processes. Individual ages are generally consistent with the accepted 
13.7Gyr ACDM age of the universe, but are sometimes in tension with it,^^ 
given our uncertainties in knowing how early on the first stars could form. 

It is interesting to note that the FB model not only adds a billion years 
to the age of the universe as measured in a galaxy, but that the percentage 
difference is larger at earlier times. For concordance ACDM r = 0.85 Gyr 
at ^ = 6.4 when distant quasars are seen, and r = 0.365 Gyr at reionization 
at z = II. By contrast for the best fit FB model t = 1.14 Gyr at z = 6.4 
and T = 0.563 Gyr at z = 11, making the universe 34% and 54% older than 
concordance ACDM at the respective epochs. 

The fact that the expansion age of the universe determined in a galaxy is 
greater than that of concordance ACDM may alleviate but not completely 
solve various aspects of the age problem. However, the fact that the age 
of the universe is 18.6Gyr at the volume average is a basic indicator that 
the whole issue of structure formation, including calibrations, needs to be 
re-examined from first principles. 
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8.5. The Hubble bubble 

Recent analysis of Sncia data by Jha, Ricss and Kirshncr'"''* confirms an 
effect which has been known about for some time, and has been interpreted 
as our Hving in a local void,^^ the "Hubble bubble" . If one excludes Snela 
within the Hubble bubble at redshifts z < 0.025, then the value of the 
Hubble constant obtained is lower by 6.5 ± 1.8%.'^^ It is for this reason 
that Snela at z < 0.023 have been excluded from the RiessO? gold data 
set.*^ The Hubble bubble is very problematic for the standard dark energy 
paradigm, as it can contribute a systematic error to the estimate of the 
equation of state parameter, w, which is much larger than the statistical 
uncertainties that precision cosmologists hope to attain. ^^'^^ 

In the FB model the Hubble bubble is expected as a feature that re- 
sults from the implicit resolution of the Sandage-de Vaucouleurs paradox. 
Since the bare Hubble parameter characterizes the uniform "quasilocally 
measured" Hubble flow, eq. (17) also quantifies the apparent variance in 
the Hubble flow below the scale of homogeneity. The present bare Hubble 
constant, H^, is lower than the global average, H^. It represents the value 
we as observers in galaxies would obtain for measurements averaged solely 
within the plane of an ideal local bubble wall, on scales dominated by finite 
infinity regions. For us in particular, measurements of the Hubble constant 
towards the Virgo cluster would represent a scale over which we might hope 
to detect the bare Hubble constant, ff~48km sec~^Mpc~^, as this ap- 
pears to be the direction that most closely represents a local bubble wall. 
Ultralocal minivoids^^ complicate any empirical measurement. 

'Local' measurements across single voids of the dominant 30h~^ Mpc 
diameter, probe the scale over which photons on null geodesies encounter 
the fewest finite infinity domains. Such measurements should give a Hubble 
'constant' which exceeds the global average iJ^ by an amount commensurate 
to ^ voids are dominant by volume, an isotropic average will 

generally produce a Hubble 'constant' greater than for local averages 
until we sample on large enough scales that the volume average of walls 
and voids is the same as the global one. Thus the average will steadily 
decrease from its maximum at ^ 30/i~^ Mpc until the scale of homogeneity 
(~ 10Qh~^ Mpc) is reached: the Hubble bubble feature. This expectation 
is dramatically confirmed by the data points in Fig. 5, reproduced from a 
recent paper of Li and Schwarz.^^ With h = 0.617, the maximum Hubble 
'constant', up to about 20% higher than average, shoifld be attained at a 
scale of about 48Mpc, thereafter steadily decreasing until leveling out at 
the scale of homogeneity, which must be reached before the BAO scale of 
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about 168 Mpc. 

Fig. 5 also illustrates how the resolution of the Sandage-de Vaucouleurs 
paradox gives unique predictions. The two thick lines in Fig. 5 show Li and 
Schwarz's estimates of the bounds between which the data should lie us- 
ing Buchert averaging, but not accounting for the clock effect. Given that 
the data indicates a consistently higher Hubble constant below the scale of 
homogeneity, the only explanation for the lack of great statistical scatter 
between the two lines - if all clocks are effectively synchronous - is that by 
a statistical fluke we happen to be in a large local void which is expand- 
ing faster into the surrounding medium. ^^'^^ However, if the conventional 
FLRW assumptions apply beyond this Hubble bubble, then by the standard 
structure formation scenarios such a circumstance seems impossible. 

In the FB model the Hubble bubble is not a statistical fluke, but like 
cosmic acceleration is an apparent effect that arises from clock rate variance. 
Distant observers in galaxies beyond our own Hubble bubble will also detect 
a Hubble bubble centred on their location. The theoretical derivation of 
the curve that should pass though the data points in Fig. 5 will provide an 
important test of the FB model. To determine this curve, and its variance, 
requires some knowledge of void statistics, and probably some Monte-Carlo 
simulations of the manner in which the scale of homogeneity is filled on 
average. Particular anisotropics which look like "large" voids^^ might be 
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expected from the variance in alignments of dominant 30h~^ diameter voids 
below the scale of homogeneity. A systematic error of 1-2% in the CMB 
dipole subtraction^^ may well be the consequence of a Rees-Sciama dipole 
resulting from such foreground inhomogeneities.^^ 

While the decades long debate among astronomers about the value of 
the Hubble constant has been mainly dominated by arguments about sys- 
tematics,^^'^® the question of the scale of averaging has no doubt also played 
a role. Lower values of will be expected if one specifically selects direc- 
tions and scales which approximate our own bubble wall.^" However, dis- 
tance determinations on scales < 50 Mpc will generally give higher values 
of Hq. Since many of the first steps on the cosmic distance ladder are cali- 
brated on such scales, many related issues require careful reconsideration. 

8.6. Prospects for future cosmological tests 

The FB model provides a strong new competitor to the standard ACDM 
cosmology, as Table 2 shows. Indeed, since it steps out of a paradigm in 
which everything rests on a single equation of state, it will provide new and 
unique predictions. Current observational programs of course focus on the 
measurement of w{z). The difference in expectations of the FB model for 
those programs will be highlighted in a forthcoming paper. ^ 



Table 2. Model comparison 



Observation 


ACDM 


Fractal bubble model 


Snela luminosity distances 


Yes 


Yes 


BAO scale (clustering) 


Yes 


Yes 


Sound horizon scale (CMB) 


Yes 


Yes 


Doppler peak fine structure 


Yes 


Still to calculate 


Integrated Sachs-Wolfe effect 


Yes 


Still to calculate 


Primordial ^Li abundances 


No 


Yes 


CMB eUipticity 


No 


[Yes] 


CMB low multipole anomalies 


No 


[Foreground void(s): 






Rees-Sciama dipole] 


CMB Sunyaev-Zel'dovich signal83,64 


No 


Still to calculate 


Hubble bubble 


No 


Yes 


Nuclcochronology of globular clusters 


Tension 


Yes 


X-ray cluster abundances 


Marginal 


Yes 


Emptiness of voids 


No 


[Yes] 


Sandage-de Vaucouleurs paradox 


No 


Yes 


Coincidence problem 


No 


Yes 



Note: * Square brackets indicate cases where the present indication is sug- 
gestive, but much detailed calculation remains to be done. 
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The quantities that remain to be calculated in Table 2 relate largely to 
the CMB. Detailed construction of the Doppler peaks, to enable compar- 
ison to WMAP, and determination of parameters such as cTg, is a matter 
of urgency. It is a rather non-trivial exercise in recalibration of standard 
quantities, in which all steps need to be carefully reconsidered. Detailed 
quantitative calculations of CMB ellipticity and the integrated Sachs- Wolfe 
effect can only be performed in conjunction with such an analysis. 

While the calculation of the integrated Sachs Wolfe (ISW) effect will 
differ in the FB model, one must be careful not to base intuition on that of 
perturbations on FLRW backgrounds. In particular, the observed ISW sig- 
nal is believed to confirm dark energy since its consequence is a suppression 
of the gravitational collapse of matter at relatively recent times. ''"'^ If one 
replaces the words "dark energy" by "voids" in the standard qualitative 
explanation, then a probable description of the ISW effect in the FB model 
emerges. The same correlation of the ISW signal with clumped structures*'''' 
is expected; what is important is the magnitude of the effect. Mattsson^^ 
has recently proposed an extension of the Dyer-Roeder formalism, which 
may have some relevance for such calculations. 

In offering a new paradigm for cosmology, every observational test nat- 
urally has to be revisited from first principles. Thus there are potential 
tests which relate to both strong and weak gravitational lensing, for ex- 
ample, which have not been listed in Table 2. Qualitatively one expects 
voids to be emptier in the FB model than in LCDM structure formation 
simulations. However, the specification of the background average surfaces 
of homogeneity, in terms of a uniform bare Hubble expansion gauge, and a 
correct post-Newtonian approximation on such a background, need to be 
examined before structure formation simulations are attempted. 

Ultimately, apparent variance in the Hubble flow below the scale of ho- 
mogeneity will give predictions which cannot be reproduced in the FLRW 
scenario. Not only should we determine the average curve that passes 
through the data points in Fig. 5, we should ultimately collect thousands 
of Snela measurements or other distance measurements on scales up to 
200Mpc, and test the correlation of the apparent Hubble flow with actual 
structure: a sort of Hubble flow tomography of the nearby universe. 

9. Conclusion 

A true "concordance cosmology" should agree with all reliable observations, 

and not just a carefully selected subset. A glance at Table 2 reveals that 
there are many anomalies in the standard ACDM model. Much attention 
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has been focused on the CMB anisotropics, on account of the spectacular 
success of the WMAP mission. However, it must be recalled that many 
ingredients go into the analysis of the CMB temperature fluctuations. In 
looking at these ingredients we must ask "what is the weakest link"! 

In my view, the weakest link is not primordial nucleosynthesis,^^ which 
is based on nuclear physics that we understand very well, but the cos- 
mological model. The weakest link comes from abandoning the theoretical 
principles that wc arc careful to apply in other circumstances. In partic- 
ular, in general relativity one should model the universe with the matter 
distribution one observes, rather than trying to impose onto the universe a 
simple mathematical solution based on some simplifying symmetry. 

It is unfortunate that general relativists have been obsessed by ex- 
act solutions of Einstein's equations, whether they involve likely or un- 
likely approximations for the matter distribution. We should face up to the 
fact that the solution for the actual matter distribution is analytically in- 
tractable, and therefore the question of cosmological averaging'''^'^^'^^'^*"'^° 
is paramount. Furthermore, once we do take an average we must address the 
fundamental problem that the relationship of rods and clocks at one point 
to those at a distant point, a conceptual centrepiece of general relativity, 
is highly non-trivial once gradients in spatial curvature and gravitational 
energy are considered. In an expanding universe these involve subtle dy- 
namical aspects of general relativity, which cannot be localized at a point 
on account of the equivalence principle. 

Some colleagues when presented with Table 2, and the fact that the 
right hand column is based entirely on general relativity with no new or 
exotic physics (beyond a need for non-baryonic dark matter), suggest I 
should invoke Ockham's Razor at this point and declare the cosmological 
constant dead. Other colleagues, who sometimes confess to Newtonian in- 
tuition, are of the view that a clock-rate difference of 38% accumulated 
between bound systems and the volume average over the lifetime of the 
universe is nonetheless so great that I must be mistaken somewhere, in 
spite of Table 2. 

In my view caution should always be exercised, but this includes caution 
with the conceptual basis of our theory and the operational interpretation 
of measurements. To those who are uncomfortable with my proposal about 
cosmological quasilocal gravitational energy let me ask the following: With- 
out reference to an asymptotically flat static reference scale, which does not 
exist given the universe is expanding, and without reference to a background 
which evolves by the Priedmann equation at some level, an assumption 
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which is manifestly violated by the observed inhoniogeneities, What keeps 
clocks synchronized in cosmic evolution? Please explain. 

In the FB proposal the fact that our FLRW approximation has served so 
well is understood as a consequence of the fact that the bare Hubble flow is 
uniform, despite large-scale inhomogeneitics. Uniformity of the quasilocally 
measured Hubble flow does not imply uniformity of spatial curvature, nor 
of gravitational energy, however. If we examine the Hubble flow over scales 
on which the average gradients arc not large, we will not sec large statistical 
scatter in the Hubble flow. But averaged over larger scales, below the scale of 
homogeneity, we will see a variance in the apparent Hubble flow, that agrees 
with observations of the statistical properties of voids, "'^^ and which seems 
to accord with Fig. 5. The fact that there is a statistically dominant void 
scale^^ may well be associated with the evolution of the scale corresponding 
to the second (rarefaction) Doppler peak in the CMB anisotropy spectrum, 
just as evolution of the first (compression) peak appears to statistically 
define the cutoff for the scale of homogeneity. 

Any proposal which seeks to genuinely shift the paradigm against ac- 
cepted intuition, even a conservative proposal based on a principled ap- 
proach^'* to our best theory of gravitation, naturally faces much initial 
scepticism. Nonetheless, I suspect that in future cosmologists will wonder 
how we ever could have expected that two decades of precision measure- 
ments would simply confirm the naive cosmological models constructed in 
the 1920s, long before we knew what the universe actually looked like. I 
suspect that the coming decades will see us explore and finally understand 
the arcane territory of general relativity relating to quasilocal gravitational 
energy,^^'^^'^^'^^ informed directly by new cosmological observations. 
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